

// Plot the distribution of RS estimates 

use $rspath/rs_hs02_19962006.dta, clear 

g top=rs_hs6_b+1.64*rs_hs6_se
g bottom=rs_hs6_b- 1.64*rs_hs6_se
sum  rs_hs6_se rs_hs6_b, d
g out=(bottom<0 & bottom!=.) | (top!=. & top>5)
sort rank
gr tw (rspike top bottom rank if out==0, ls(ci) lw(vvthin) lc(gs14))  (line rs_hs6_b rank, lp(thick) lc(gs1)), legend(off) xtitle("Rank - relationship stickiness") ytitle("Relationship stickiness")
graph export $outputpath/figure1.pdf, replace



